Decline of the space density of quasars between z = 2 and z = 4 
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ABSTRACT 

We define a new complete sample of 13 optically-luminous radio quasars (M J 4s(1450A)< 
-26™9 and logPi^GH^WHzr 1 ) > 25.7) with redshift 3.8 < z < 4.5, obtained by cross- 
correlating the FIRST radio survey and the APM catalogue of POSS-I. We measure the space 
density to be (1.0 ± 0.3) Gpc -3 , a factor 1.9±0.7 smaller than the space density of similar quasars 
at z fv 2 (FBQS). Using a new measurement of the radio-loud fraction of quasars we find that at 
z = 4 the total space density of quasars with M A s(1450A)< -26™9 is (7.4 ±2.6) Gpc" 3 . This is 
a factor 1.8 ± 0.8 less than the space density at z fa 2, found by the 2dF quasar survey. This z = 
2 / z = 4 ratio, consistent with that of the radio-loud quasars, is significantly different from the 
ratio ~ 10 found for samples including lower-luminosity quasars. This suggests that the decline 
of the space density beyond z « 2 is slower for optically-luminous quasars than for less-luminous 
ones. 

Subject headings: cosmology: observations, early-universe, galaxies: evolution, galaxies: formation, 
quasars: general, surveys ^ 



1. Introduction 

The space density of bright quasars is a factor 
« 100 higher at z — 2 than at z = 0. Boyle et 
al. (2000; the 2dF survey) have shown that this 
increase can be described in terms of pure luminos- 
ity evolution, i.e. the shape of the luminosity func- 
tion does not change with z, up to at least z = 2.3. 
At higher redshifts evidence for a decrease in space 
density, consistent with pure luminosity evolution, 
has been presented by several authors (e.g. Ken- 
ncfick et al., 1995 and references therein, but see 
also Irwin et al. 1991). The Sloan Digital Sky 
Survey (SDSS, Fan et al. 2001a) has pushed ex- 
ploration of the quasar phenomenon out to z > 6, 
the epoch of formation of the first massive black 
holes (and galaxies) (Becker et al. 2001b). The 
space density of quasars at z ~ 4 with M J 4s(1450 
A) < -25 m 5 (n = 1 and A = 0, H = 50 km s" 1 
Mpc" 1 , as used throughout this paper) was found 
by Fan et al. (2001b) to be more than one order 
of magnitude lower than that found by Boyle et 
al. (2000) at z - 2.3. Fan et al. (2001b) also 
found substantial deviations from pure luminosity 
evolution: the bright end of the quasar luminos- 
ity function at z > 4 appears to be flatter than 
found at lower redshifts, i.e. the decrease in space 
density with increasing redshift depends on lumi- 
nosity. 

Radio-loud quasars are a small subset of the 
quasar population, but as stressed by several au- 
thors (e.g. Peacock 1985; Dunlop & Peacock 1990; 
Shaver et al. 1996, 1999), selection in the radio 
drastically reduces contamination by foreground 
stars with optical colors similar to high-redshift 
quasars, and is unaffected by both intrinsic and ex- 
trinsic dust obscuration and reddening. Thus the 
high-redshift behaviour of the radio-loud quasar 
luminosity function is a robust probe of the abun- 
dance of massive objects at very early cosmolog- 
ical epochs. A second important motivation for 
radio searches of high-z quasars is understanding 
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the origin of the radio phenomenon. Differences, 
if any, in the early evolutionary behaviour of radio 
and optical activity would provide valuable clues. 

Several searches for high-redshift radio-loud 
quasars have been carried out (e.g.: Shaver et 
al. 1996, 1999; Hook et al. 1998a; Jarvis & 
Rawlings 2000; Jarvis et al. 2001; Snellen et al. 
2001; Hook et al., 2002). Shaver et al. (1996, 
1999) obtained complete redshift information (i.e. 
no optical magnitude limit) for a sample of 442 
radio-loud flat-spectrum quasars with flux densi- 
ties S2.7GHZ > 0.25 Jy. They argued for a drop 
by more than a factor of 10 in space density from 
z ~ 2.5 to z ~ 6, consistent with that found in 
the optical. 

Hook et al. (1998b) reported indications that 
the peak in the quasar space density may have 
occured at earlier cosmological epochs for more 
radio-luminous sources, the decline in space den- 
sity with redshift being less pronounced for these 
radio sources. Jarvis & Rawlings (2000) ana- 
lyzed the top decade in radio luminosity (log 
-P2.7GH«(WHz _1 ) > 28, quasars and galaxies) 
of sources from the Parkes Half-Jansky Flat- 
Spectrum Sample. They deduce a decline by a 
factor ~ 4 between z ~ 2.5 and z ~ 5, less 
abrupt than envisaged by Shaver et al. (1996, 
1999) and consistent with that measured by Dun- 
lop & Peacock (1990) for radio sources in gen- 
eral. Jarvis et al. (2001) investigated the high- 
redshift evolution of the most-radio-luminous 
(28.8 < logPisiMH^WHz" 1 ) < 29.8) steep- 
spectrum radio sources using three samples se- 
lected at low frequencies (3CRR, 6CE, and 6C*). 
They found the data consistent with a constant 
space density between z ~ 2.5 and z ~ 4.5 and 
exclude a decline as steep as suggested by Shaver 
et al. (1996, 1999). These results are not neces- 
sarily inconsistent, since different ranges of radio 
luminosity and different frequency selections were 
used. 

No evidence for a rapid decline of the space den- 
sity of soft-X-ray selected quasars at z > 2.7 was 
apparent in the ROSAT survey data analyzed by 
Miyaji et al. (2000). 

In this paper we report the results of a search 
for z i=a 4 radio-loud quasars, down to radio 
flux densities thirty times fainter than in previous 
searches, and with no selection on radio spectral 
index. The paper is organized as follows: 
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In Sect. 2 we describe the construction of our 
sample of z = 4 quasars. 

In Sect. 3 we calculate the space density at z w 4, 
and compare it with that at z w 2 derived from 
a sample of quasars of similar optical and radio 
luminosity. 

In Sect. 4 we calculate the z w 4 space density of 

all quasars optically brighter than M j 4b(1450A) < 

— 26™9, using a new determination of the fraction 

of radio loud quasars at that redshift, and compare 

it with that measured at z = 2. 

Section 5 contains the discussion, and Sect. 6 our 

conclusions. 

2. Sample of z w 4 radio quasars 
2.1. Candidate selection 

The broad-band optical color of quasars red- 
dens dramatically with redshift, as the Lyman-a 
forest and the 912-A Lyman-limit leave the blue 
band (Fig. 1). To obtain a sample of high-redshift 
quasars, we therefore sought optical identifications 
of FIRST radio sources (Becker et al. 1995, White 
ct al. 1997) with objects catalogued by the Auto- 
mated Plate Measuring facility (APM; (McMahon 
& Irwin 1992) as being star-like on the POSS-I E 
(red) plates, with red magnitude E < 18™8, and 
with measured color O — E > 3, or O fainter than 
the plate limit (O ~ 21 m 7). 

The FIRST radio catalogue (version of July 
2000) includes w 660,000 radio sources with peak 
flux density Siaghz > 1 nijy in 2.24 sr and 
7 h < RA < 17 h , -5° < DEC < 57°, of which 
2.138 sr are covered by the APM catalogue. 

We sought optical identifications lying closer 
than 1'.'5 to each radio position, corresponding to 
~ 3 standard deviations in the expected combined 
radio and optical positional errors (see Fig. 2). 

Searching to a larger radius would have netted 
1% more identifications, but the contamination 
from foreground stars rises rapidly with search ra- 
dius (see again Fig. 2). A similar search radius 
(T/2) was used by White et al (2000) in their 
search for First Bright Quasar Survey (FBQS) 
candidates. Our search yielded a starting list of 
2323 candidates. 

Trial spectroscopy showed that w 85% of these 
'stellar' candidates are faint galaxies. 

In order to reduce contamination by these 




Fig. 1.— POSS-I O-E color vs redshift for high- 
redshift quasars (a) with 3.6 < z < 4.5, E < 18™8, 
from Djorgovski (2001, z > 4) and from NED (as of 
May 2002, z < 4); (b) the 18 high-z FIRST quasars 
identified in our search (Section 2.2). Filled sym- 
bols indicate measured O — E colors, arrows are 
lower limits (e.g. > 3 for E < 18 m 7, O limit 
21™7). The quasars in (a) were selected in part 
on the basis of color, but nevertheless trace the 
typical reddening of a quasar with increasing z 
as the continuum redward of the Lya line moves 
out of the observed O band (red limit w 5000 A; 
prediction for typical quasar shown as solid curve, 
computed in a similar way as described in the Ap- 
pendix). As redshift increases through z w 4, O- 
band shifts below the 912 A Lyman-limit and the 
O — E color rises rapidly. As redshift increases 
above 4.5, the predicted brightness of a quasar in 
E band drops sharply, so few are likely to be de- 
tected. The dashed line indicates z — 3.8, above 
which we believe selection by color to be nearly 
complete (Sect. 2.3). 
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galaxies, we ran the SExtractor program (Bertin 
& Arnouts 1996) on images of the 2323 candidates 
on the POSS-II red plates, which have better sensi- 
tivity and spatial resolution than POSS-I. Eighty- 
two percent of the candidates (1901 objects) were 
clearly extended (i.e. galaxies). Seventy-three 
candidates could not be found by SExtractor, ei- 
ther because they were defects on the POSS-I 
plate (62 objects) or due to confusion effects near 
a bright star (11 objects). 

The remaining 349 candidates were sought in 
the Minnesota APS catalogue of POSS-I (Pen- 
nington ct al. 1993). Fifty-five of the candidates, 
undetected by APM on the O plate (i.e. poten- 
tially very red), were detected by APS. These have 
O — E < 3, so were removed from the sample. 
This left 294 candidates, 194 of which we classi- 
fied, on the basis of the SExtractor-generated pa- 
rameters, as 'star-like', and 100 as 'possibly star- 
like'. CCD images of a third half of these were 
obtained in l'.'3-seeing with the Loiano 1.5-m tele- 
scope of the Astronomical Observatory of Bologna. 
On these images, 25 of the 70 'star-like', and 28 
of the 33 'possibly star-like' candidates were ex- 
tended. We have not investigated further the 100 
'possibly star-like' objects. As a result we expect 
an incompleteness of m 100 x 5/33 = 15(±4) stel- 
lar candidates. 

2.2. Spectroscopy 

Spectra were obtained of 121 of the 194 can- 
didates, using the IDS spectrograph of the Isaac 
Newton Telescope (Benn ct al. 2001, Holt et al., 
in prep.). Eighteen of the 121 (15%) are quasars 
with 3.6 < z < 4.5. The fraction of high-redshift 
quasars varies little with the optical magnitude 
of the candidates, but depends strongly on radio 
flux density: (7.7 ± 2.9)% for S 1A gh z < 10 mJy 
and (37 ± 11)% for Si.4g.Hz > 10 mJy. Amongst 
the 73 candidates not observed spectroscopically, 
there are 6 with Siaghz > 10 mJy and 67 with 
Siaghz < 10 mJy. On the basis of the above high- 
rcdshift-quasar fractions, and accounting also for 
the ~ 15 missed stellar candidates (Sect. 2.1), we 
expect to have missed a total of 9 ± 2 high redshift 
quasars, i.e. the search is « (67 ±5)% complete 
or, equivalcntly, corresponds to an area on the sky 
of 1.432 sr. 
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Fig. 2. — The distributions of radio-optical separa- 
tions r for the star-like candidates classified spec- 
troscopically as quasars, galaxies and stars (total 
194 objects). For quasars and galaxies, the dashed 
curves are of the form r • exp{— r 2 /2er 2 ), with a = 
Of.' 5. For stars, the dashed straight line indicates 
the form expected for misidentifications with fore- 
ground stars. 
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2.3. Sample completeness 

Our search for high-redshift quasars within this 
area will be incomplete for a number of other rea- 
sons, that we discuss below: 

1) Incompleteness of the APM catalogue. Ac- 
cording to McMahon et al. (2001), APM classifies 
correctly 92% of E < 19 m objects on at least one 
of the O or E plates. Assuming that the APM E 
and O plate misclassifications are independent and 
occur at a similar rate, then a 92%-correct classi- 
fication on at least one plate implies 72%-correct 
classification on individual plates. This last figure 
is a lower limit on the incompleteness of our E 
plate based star-like classification (the classifica- 
tion on the O and E plate cannot be completely 
independent). 

For an independent estimate of our incomplete- 
ness, we searched for APM detections of 100 z > 4 
quasars with E < 19 m from the list of Djorgovski 
(2001), which were originally found either on the 
POSS-II plates or in SDSS (CCD based). All 100 
were detected by APM. Eighteen of these quasars 
appeared extended in the APM E band, i.e. the 
search was (82 ± 4)% complete. 

As a further check we searched the APM cata- 
logue for counterparts of 4949 objects with 16 m < 
E < 19 m catalogued as star-like in the APS. APM 
classifies as extended (13 ± 0.5)% of these objects 
and misses a further (3 ± 0.3)%, i.e. is 84% com- 
plete. Because the APS classification algorithm 
(neural network) is believed to be more reliable 
than that of APM, we assume that disagreements 
about classification by APM and APS are due 
mainly to errors in the former. 

To check the APS classification reliability, we 
constructed a sample of 190 quasars, 140 from the 
FBQS survey with 16 m < E < 18 m and 50 from 
the B3VLA QSS survey (Vigotti et al. 1997) with 
18 m < E < 19 m . Only 10 quasars were classified 
extended by APS, i.e. an incompleteness of 5%. 
This is three times smaller than the APM incom- 
pleteness, consistent with our assumption above. 

The above two estimates of APM completeness 
are in good agreement and we assume a value of 
(84 ± 2) %. In addition, a few star-like objects 
might have been misclassified extended when we 
analysed the images with SExtractor, or at the 
stage of CCD imaging, but the number is likely to 



be very small. 

2) Omission of quasars with extended radio 
emission. Optical counterparts were sought at the 
positions of individual FIRST radio sources, but 
not at the mid-points of possible double sources. 
In the First Bright Quasar Survey (FBQS) in the 
South Galactic Cap, Becker et al. (2001a) also 
sought counterparts along the lines joining possi- 
ble double sources of separation up to 30", and 
found that a search at the positions of individual 
radio components is 96% complete. At z ~ 4, the 
rest-frame frequency of observation is even higher 
than for FBQS, favouring more flat-spectrum com- 
pact sources. We therefore expect our sample to 
be more complete, and we assume a completeness 
of (98 ± 1)%. 

3) Incompleteness due to errors in the radio- 
optical positions. The distribution of the radio - 
optical separations of the quasars (Fig. 2) is con- 
sistent with a rms error of Of! 5, implying that our 
search out to radius l'/5 is ~ 99% complete. 

4) Incompleteness due to color selection. The 
fraction of z > 3.8 quasars with O — E > 3 m can 
be estimated approximately from the color distri- 
bution in Fig. la (taken from the literature). 81 of 
the 85 quasars with z > 3.8 have O — E > 3, i.e. 
this color criterion finds (95 ±1.5) % of quasars 
with z> 3.8, but the efficiency falls rapidly with 
decreasing redshift. We therefore restrict our sam- 
ple to z > 3.8. 

Combining the above, the completeness of our 
search for z > 3.8 quasars in FIRST is « (77±6)% 
over the effective area (Sect. 2.2) of 1.432 sr. 

5) Incompleteness of the FIRST radio survey. 
The FIRST catalog is severely incomplete at faint 
flux densities. From Prandoni et al. (2001) we 
estimate a completeness of « 50% for £ < 1.25 
mJy and of w 75% for 1.25 mjy < S < 1.5 mjy 
This flux-density-dependent incompleteness is in 
addition to that discussed above. 

2.4. The z s=a 4 radio-loud sample 

The 13 radio-loud quasars with z > 3.8, the 
QziRL sample, are listed in Tab. 1. 

The E magnitudes were recalibrated using APS 
catalogue (rms error 0™2, better than APM, be- 
cause the APS photometry was calibrated on a 
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plate-by-plate basis) and corrected for Galactic ex- 
tinction. After re-calibration the limiting magni- 
tude of the sample is E = 19 m 

The epoch-1950 absolute magnitudes in the 
rest-frame UV continuum, Mab (1450 A), were ob- 
tained from the E magnitudes, using a template 
quasar spectrum (the procedure is described in the 
Appendix). The 13 quasars are optically very lu- 
minous, Mab (1450 A)< -26 m 9. 

The radio luminosities were calculated using 
a median spectral index cur = —0.3, rms = 0.4 
(S„ oc v a ), and range logPi. 4G Hz(W Hz" 1 ) > 25.7. 
The median spectral index was obtained via 5- 
GHz observations of some sources using the Effcls- 
berg radio telescope (Holt et al. 2003, in prep.). 

As a check of our method of deriving Mab (1450 
A), we made, for 8 of the quasars, independent 
measurements based on recent /-band photome- 
try obtained during service time at the 1.0-m Ja- 
cobus Kaptcyn Telescope in La Palma (Holt et 
al., in prep.). At z = 4.2, I band corresponds 
to rest-frame wavelength <~ 1540A, close enough 
to the reference wavelength 1450A that we can 
K-correct using the median optical spectral index 
o-opt = -0.7 from Fan et al. (2001c). For two of 
the 8 quasars, 0831+52 and 1423+22, we found a 
large difference between the Mab values derived 
using these two methods, the /-photometry giving 
Mab respectively dimmer by 0™4 and brighter by 
1™0. These differences arc probably due to intrin- 
sic changes in luminosity between epochs 1950 and 
2002. For the remaining 6 quasars the mean differ- 
ence is 0.0 ± 0.1, suggesting that any systematic 
error in our Mab magnitude scale is < 0™1. 

The mean radio and optical luminosities of 
the sample are: (logPnGHz) = 26.6 ± 0.2 and 
(Mab) = -27^6 ± m 2, the median values being 
26.8 and -27^5 respectively 

3. Decline of the space density of radio 
loud quasars between z = 2 and z = 4 

3.1. Comparison sample of z ~ 2 radio 
quasars 

For a sample of radio quasars with comparable 
radio and optical luminosities, but at z m 2, we 
draw on the 0.817-sr First Bright Quasar Survey 
(FBQS, White et al. 2000). FBQS uses O, E mag- 
nitudes recalibrated using APS, with a limiting 



magnitude E — 17^8. Of the original 1238 quasar 
candidates, 636 were confirmed quasars. No spec- 
trum is available for 108 of the 1238 candidates, 
i.e. the spectroscopy is 91% complete, correspond- 
ing to an effective area on the sky of 0.743 sr. 

Of the 45 FBQS quasars with 1.8 < z < 2.4, 
23 have luminosities in the range of the Q z arl 
sample, i.e. log Pi. 4 GHz > 25.7 and Mas (1450 
A)< -26 m 9. 

The mean radio and optical luminosities val- 
ues are (logPi. 4G ^) = 26.6 + 0.2 and (Mab) = 
— 27^5 ± 0.1. The median values are 26.4 and 
— 27^4 respectively (consistent, within the statis- 
tical noise, with those of the z = 4 sample). We 
thus have a lower-redshift sample directly compa- 
rable with our high-redshift one. Fig. 3 shows 
the distributions of the sources in radio and opti- 
cal luminosity and redshift for each sample, for 
the Einstein - de Sitter cosmology and for the 
concordance model (see Sect. 3.2) 

This sample (Tab. 2) will be referred as the 
Qz2RL sample. The radio luminosities Piaghz 
were calculated assuming a radio spectral index 
cir = —0.3. This is the median spectral in- 
dex obtained for a subsample with GB6 4.8-GHz 
flux densities (Gregory et al. 1996), and is sim- 
ilar to that of the z ~ 4 sample. The absolute 
magnitudes Mas (1450 A) have been calculated 
from the E magnitudes using a K-correction = 
2.5 log(l + z )( 1 +"opt) and aQpT = _ .5, that of 

the FBQS composite spectrum (Brothcrton et al. 
2001). 

As an independent check, we also calculated the 
absolute magnitudes using the quasar template 
(following a procedure similar to that used at z = 
4, see Appendix), obtaining similar results, within 
±0 m l. 

This sample will suffer from some of the same 
incompletenesses as the Q z >±rl sample (Section 
2.3). The search for star-like objects E < 17 m 8 
in the APM will be w (96 ± 1)% complete. An 
additional 4% of objects will be missed due to 
not searching at the mid-points of possible dou- 
ble radio sources (Becker et al. 2001a). No losses 
are expected due to color selection, nor for in- 
completeness of the FIRST at faint flux densities 
(these sources are brighter than 3 mJy). In to- 
tal, the Qzirl sample of quasars is expected to be 
(92 ± 1.4)% complete over 0.743 sr. 
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Table 1: Properties of the z = 4 radio- loud sample 





RA 




Dec 


E 
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- E 


z 


SlAGHz 


\ogP 1AGHz 


M AB 


Zmax 


P 






J2000 










mJy 


WHz -1 


(14501) 




Gpcr 3 


07 


25 18 


.29 


+37 05 17.9 


18.7 


> 


3.4 


4.330 


26.6 


27.11 


-26.9 


4.379 


0.0517 


07 


47 11 


.16 


+27 39 03.6 


17.8 


> 


4.5 


4.170 


1.1 


25.70 


-27.6 


4.366 


0.0530 


08 


31 41 


.68 


+52 45 17.1 


15.8 




3.9 


3.870 


1.3 


25.72 


-29.8 


4.500 


0.0432 


08 


39 46 


.11 


+51 12 03.0 


18.9 


> 


3.2 


4.420 


41.6 


27.32 


-27.3 


4.438 


0.0472 


09 


18 24 


.38 


+06 36 53.5 


19.0 


> 


3.0 


4.180 


26.5 


27.09 


-26.9 


4.206 


0.0731 


09 


41 19 


.36 


+51 19 32.3 


18.0 


> 


4.1 


3.850 


2.5 


26.00 


-27.6 


4.470 


0.0450 


10 


53 20 


.42 


-00 16 50.4 


17.8 


> 


3.6 


4.289 


13.8 


26.82 


-27.7 


4.483 


0.0441 


10 


57 56 


.27 


+45 55 53.0 


16.9 




3.9 


4.110 


1.1 


25.70 


-28.5 


4.350 


0.0544 


12 


11 34 


.41 


+32 26 15.8 


18.2 


> 


3.0 


4.097 


3.7 


26.21 


-27.2 


4.421 


0.0484 


13 


09 40 


.61 


+57 33 09.1 


18.3 


> 


3.6 


4.274 


11.3 


26.73 


-27.1 


4.406 


0.0495 


13 


25 12 


.52 


+ 11 23 29.9 


18.6 


> 


2.9 


4.400 


71.1 


27.55 


-27.5 


4.455 


0.0460 


14 


23 08 


.19 


+22 41 57.5 


18.7 


> 


3.2 


4.316 


35.4 


27.23 


-26.9 


4.377 


0.0519 


16 


39 50 


.51 


+43 40 03.6 


17.7 




3.9 


3.990 


25.2 


27.03 


-27.8 


4.500 


0.0432 



Table 1: The columns give: (1 - 2) J2000 optical coordinates, (3) E magnitudes recalibrated using APS; 
(4) APM O — E color; (5) spectroscopic redshift; (6) FIRST radio flux density (peak or integrated, see 
next column), (7) radio luminosity, calculated using the peak flux density for Siaghz < 10 mJy and the 
integrated flux density for Siaghz > 10 mJy, as recommended by Ivezic et al. (2002, see their Fig. 4), and 
K-corrected using the median radio spectral index an = —0.3 (S„ <~ v a , as used throughout this paper), 
(8) absolute magnitude Mab at rest-frame 1450 A (see Appendix), (9) maximum redshift out to which the 
quasar could have been detected (i.e. with E < 19™0, Siaghz > 1 mJy), (10) contribution of this source to 
space density, assuming effective area 1.432 sr, but without correction for the incompletenesses discussed in 
Sect. 2.3 (see text). 



3.2. The space density of radio-loud quasars 
at z m 4 and at z w 2 

We have two shells of volume: one defined by 
1.8 < z < 2.4 (which has been chosen to yield 
a good statistical sample around z — 2) and the 
other by 3.8 < z < 4.5 (determined by the fact 
that at z larger than 4.5 the predicted brightness 
of a quasar in E band drops sharply, so few are 
likely to be detected). For each of the two sam- 
ples of quasars, the space density was calculated 
as follows. For each source we computed z max>r 
and z max ,o the maximum redshifts at which the 
source would reach the radio flux-density limit or 
the optical magnitude limit of the sample. 

In addition, for the z = 4 sample of quasars, 
we calculated the minimum redshift at which 
the quasar could have been seen, since the K- 
correction could in principle cause the quasar to 
dim with decreasing redshift (Fig. 6). For all 
the sources we found this minimum redshift to be 
< 3.8. 



The maximum redshift z max at which the 
quasar could have been included in the sample 
is the minimum of z max , r and z max , Q . 

Just two of the quasars in the z = 4 sample are 
radio limited {z max . r < z maX:0 ), while in the z = 2 
sample there are none. The uncertainty in z max , r 
due to the assumption of an average radio spectral 
index (—0.3) is ~ 0.1 and implies a negligible effect 
on the calculated space density (< 2%). 

The lower limit on redshift is the sample limit zq 
(1.8 for the z«2 sample, 3.8 for the z m 4 sample). 
The space density contributed by each quasar is 
the inverse of the volume of a shell bounded by Zq 
and z max , and with surface area equal to the effec- 
tive area covered by the sample (0.743 and 1.432 
sr for the Q Z 2RL and the Qz4RL samples respec- 
tively). This volume is: 
V = k((l - -T—l ) 3 - (1 - -ri ) 3 ) 

where k = 575 Gpc 3 per sr. 
Summing the contributions of individual sour- 
ces, and normalising by area (Tabs. 1, 2), we ob- 
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Table 2: Properties of the z = 2 radio- loud sample 





RA 






Dec 




E 


O-E 




z 


SlAGHz 


l0gPi.4Gffz 




z 


max 


P 








J2000 
















mJy 


WHz" 1 


(1450A) 






Gpc~ 3 


07 


29 


28. 


,47 


+25 


24 


51 


.99 


17.2 


0.6 


2 


.303 


6.8 


26.05 


-27.4 


2 


.400 


0.0761 


07 


58 


13 


.94 


+26 


24 


18 


.75 


17.1 


0.7 


1 


.915 


13.9 


26.22 


-27.1 


2 


.400 


0.0761 


08 


04 


13. 


.67 


+25 


16 


34 


.07 


17.4 


0.7 


2 


.302 


3.5 


25.76 


-27.2 


2 


.400 


0.0761 


08 


57 


26. 


.95 


+33 


13 


17 


.24 


17.1 


0.3 


2 


.338 


21.3 


26.56 


-27.6 


2 


.400 


0.0761 


09 


04 


44. 


.33 


+23 


33 


54 


.05 


16.9 


0.3 


2 


.244 


107.3 


27.23 


-27.7 


2 


.400 


0.0761 


09 


10 


54. 


.18 


+37 


59 


15 


.11 


17.2 


0.0 


2 


.158 


265.9 


27.59 


-27.3 


2 


.400 


0.0761 


10 


12 


11. 


.43 


+33 


09 


25 


.96 


17.3 


0.0 


2 


.260 


196.3 


27.50 


-27.3 


2 


.400 


0.0761 


10 


35 


10. 


.98 


+35 


10 


19 


.51 


17.3 


-0.2 


1 


.955 


22.9 


26.45 


-26.9 


2 


.400 


0.0761 


10 


36 


41. 


.94 


+25 


02 


36 


.63 


16.9 


0.7 


2 


.004 


118.5 


27.18 


-27.5 


2 


.400 


0.0761 


10 


54 


27. 


.15 


+25 


36 


00 


.33 


16.9 


1.4 


2 


.400 


3.0 


25.72 


-27.8 


2 


.400 


0.0761 


11 


22 


41. 


.46 


+30 


35 


34 


.88 


16.4 


0.4 


1 


.809 


10.0 


26.03 


-27.8 


2 


.400 


0.0761 


12 


25 


27. 


.38 


+22 


35 


12 


.72 


16.2 


0.2 


2 


.043 


3.9 


25.71 


-28.2 


2 


.400 


0.0761 


12 


28 


24. 


.98 


+31 


28 


37 


.74 


15.5 


0.4 


2 


.199 


323.8 


27.69 


-29.0 


2 


.400 


0.0761 


13 


24 


22. 


.53 


+24 


52 


22 


.25 


17.4 


1.5 


2 


.357 


4.9 


25.92 


-27.3 


2 


.400 


0.0761 


13 


45 


20. 


.42 


+32 


41 


12 


.59 


17.6 


0.3 


2 


.253 


14.6 


26.36 


-27.0 


2 


.400 


0.0761 


14 


15 


28. 


.47 


+37 


06 


21 


.49 


17.4 


0.9 


2 


.360 


409.7 


27.85 


-27.3 


2 


.400 


0.0761 


14 


16 


17. 


.38 


+26 


49 


0G 


.20 


17.7 


0.8 


2 


.297 


2.9 


25.68 


-26.9 


2 


.400 


0.0761 


14 


21 


08. 


.72 


+22 


41 


17 


.42 


16.8 


0.1 


2 


.181 


4.7 


25.85 


-27.8 


2 


.400 


0.0761 


14 


23 


26. 


.07 


+32 


52 


20 


.34 


16.3 


0.2 


1 


.903 


8.7 


26.01 


-27.9 


2 


.400 


0.0761 


14 


32 


43. 


.32 


+41 


03 


28 


.04 


17.4 


0.7 


1 


.963 


261.7 


27.51 


-26.9 


2 


.400 


0.0761 


1G 


03 


54. 


.15 


+30 


02 


08 


.88 


17.4 


0.7 


2 


.026 


54.2 


26.85 


-27.0 


2 


.400 


0.0761 


1G 


34 


12. 


.77 


+32 


03 


35 


.45 


16.6 


0.3 


2 


.341 


176.5 


27.48 


-28.1 


2 


.400 


0.0761 


1G 


51 


37. 


.56 


+40 


02 


18. 


.71 


16.9 


0.3 


2 


.316 


43.9 


26.86 


-27.7 


2 


.400 


0.0761 



Table 2: The columns are as for Tab. 1, except that in column: (4) O — E is the APM color using APS 
recalibration; (8) M J 4s(1450A) is calculated from the E magnitude (the O magnitude is closer in rest-frame 
wavelength, but may include Lya in the wing of the filter response); (9) the maximum redshift out to which 
the quasar could have been detected is calculated assuming E < 17^8, Siaghz > 1 mJy; (10) effective area 
0.743 sr is assumed. 



tain space densities of 0.65 + 0.18, and 1.75 + 0.36, 
quasars per Gpc 3 at z k 4.2 and z m 2 respec- 
tively. 

Correcting for completeness (77 ±5)% for the 
QziRL sample, multiplied by 60% for S < 1.5 
mJy), (92+1.4)% for the Q Z 2rl sample) the space 
densities become: 

PzArl = 0.99+ 0.28 Gpc~ 3 

Pz2RL = 1.90+ 0.47 Gpc- 3 
with ratio: 

Pz2Rl/PzArl = 1-9 ± 0.7. 

In calculating the errors, we also took into ac- 
count possible systematic shifts of 0™1 in the mag- 
nitude scale. 

Had we adopted the currently more fashionable 



cosmological concordance model with fl\ = 0.8, 
^matter = 0.2 and Ho = 70 km s _1 Mpc -1 , the 
ratio of the comoving volumes V(z — 4) jV(z = 2) 
would have been higher by a factor w 1.25. The 
derived radio and optical luminosities change by 
different amounts in the two redshift bins. E.g. 
Mab — — 26™9 in the Einstein - de Sitter cosmol- 
ogy becomes in the concordance model Mab = 
-27 m 16 at z = 2, and M AB = -2^34 at z = 
4. The same applies to the radio luminosities. To 
maintain comparable luminosities in the two sam- 
ples, eight of the z — 2 quasars have to be omitted. 
The density ratio p Z 2RL / Pz4RL then decreases by 
a factor ~ 0.81. 



8 



4. The space density of optically-bright 
quasars at z w 4 

In the previous section we measured the 
space density of optically bright (M J 4s(1450A)< 
-26™9) radio-loud (logPi. 4 > 25.7) quasars at 
z w 4. At this redshift the total space density 
of all optically-bright quasars (i.e. radio-loud + 
radio-quiet) is still poorly known, although SDSS 
should soon yield an accurate value. Here we 
obtain it from the space density of the Q z arl 
sample and from an accurate new measurement 
of the fraction of optically-bright quasars with 
log Pi. 4 > 25.7 (radio loud). 5 We indicate this 
fraction as Grl- 

4.1. The radio/optical fraction of quasars 

at z w 4 

From Djorgovski's (2001) list of z > 4 quasars 
we extracted a sample of 191 optically-selected 
quasars lying in the area of sky covered by the 
FIRST survey. For 22 of these we found a FIRST 
radio source within 3". The number of chance co- 
incidences will be negligible (the number of can- 
didates does not change if the search radius is in- 
creased to 10"), i.e. all of these are radio-loud 
quasars. 

We are interested in the luminosity range 
covered by the Q z arl sample, M,4s(1450 A)< 
— 26™9. From the reported magnitudes, which 
are APS or APM POSS-I or POSS-II R, we com- 
puted Mas (1450 A) for all 191 objects, with es- 
timated rms and systematic errors both < 0™3. 
The magnitude uncertainties are unimportant for 
our purposes, because we will see that the frac- 
tion of radio loud quasars, Grl, changes little 
over a broad range of optical luminosity. Apply- 
ing the correction factor for the incompleteness of 
the FIRST (sec Sect. 2.3) we obtain: 

G RL {M AB {U50k) < -26 m 9) = (13.3 ±4.8)% 
G RL (M ab (1450A) > -26 m 9) = (13.4 ±3.7)% 
i.e. there is no evidence that Grl depends on 
optical luminosity. The mean value is: 

G RL (-25 m > M AB (1450A) > -29 m ) = 

5 The definition of radio-loudness is rather arbitrary. We ac- 
tually define as 'radio loud' any quasar detected in FIRST, 
i.e. Si.4qhz > 1 mjy, which corresponds to log P\AGHz > 
25.7 at the median redshift of the 2 = 4 sample. Gregg et al. 
(1996) use the very similar criterion log P\AGHz > 25.5. 



(13.4 ±3.0)%. 

This is consistent with the less accurate esti- 
mate (12 ± 6) % by Stern et al. (2000). 

In the next paragraph we will use this ratio and 
the space density of the Q z arl loud sample to find 
the space density of all quasars at z ~ 4 with 
M AB < -26 m 9. ' 

4.2. The luminosity function of optically 
bright quasars 

Above we measured a radio-loud fraction 
Grl = (13.4 ± 3.0)% for quasars at z ~ 4. We 
can therefore calculate the space densities of all 
quasars at this function of M AB (U50k), 

simply by dividing the corresponding space densi- 
ties of the radio-loud quasars by Grl e.g.: 

p{M AB (U50k) < -26 m 9) = (7.4 ± 2.6) Gpc~ 3 

and 

p{M AB (U50k) < -27 m 5) = (3.8 ± 1.7) Gpc^ 3 

These numbers are consistent with the results 
of Fan et al. (2001b) at z = 4.3 (Fig. 4). 
Their maximum-likelihood luminosity function 
gives a density of 4.1± 3 \Gpc-* at M AB {U50A) < 
-26. 9 m . 

In Fig. 4 we compare our cumulative quasar 
luminosity function with that from 2dF at z w 2 
(Boyle et al. 2000). We use the actual data 
(their Fig. 4), rather than their pure-luminosity- 
evolution model fit, which tends to overestimate 
the space densities at the bright end (both arc 
shown in Fig. 4). In order to make this compar- 
ison we have converted their M B to M AB (14:50 
A) assuming M B = M AB (1450A) + 1.2a opt + 0.12 
(Kcnncfick et al. 1995) and adopting a op t — — 0.5 
(Boyle et al. 2000). In estimating the errors we 
allow also for possible systematic differences of 0.1 
in the magnitude scales. The 2dF space density at 
z w 2 is: 

p(M AB (1450A) < -26 m 9) = (13.5 ±4) Gpc" 3 

Thus the ratio between the space densities, for 
M AB (1450A)< -26^9, at z w 2 and at z w 4 
is: 

Pz2/pzi = 1.8 ±0.8 
If we extrapolate the z = 2 luminosity function to 
brighter magnitudes we get: 

p(M j4B (1450A) < -27 nl 5) w 2.0 Gpc' 3 
and the ratio is less than 1 at Af4 B (1450A) < 
-27 m 5. 
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Fig. 3. — The plots of Radio and Optical lumi- 
nosity vs z, for the Q z arl sample (left plots) and 
for the Q Z 2RL sample (right plots). Filled symbols 
are for the Einstein - de Sitter cosmology, empty 
symbols are for the concordance model (see Sect. 
3.2). 
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Fig. 4. — The cumulative luminosity function of 
optically-luminous quasars at z <~ 4.2, derived in 
Sect. 4.2 (filled dots, 2 sample error bars shown), 
compared with SDSS results for z = 4.3 (shaded 
area, ±ler confidence level). The solid curve is the 
cumulative luminosity function obtained by inte- 
grating the 2.0 < z < 2.3 luminosity function from 
the 2dF survey (Boyle et al. 2000), the dashed 
curve is their global model best fit. 



Had we used the pure-luminosity-evolution 
model the densities would have been larger by 
a factor s=a 1.7 and « 2.0 and the density ratios 
would have been m 3.0 and w 1. 

The change in space-density ratio with optical 
luminosity is evident in Fig. 4. 

5. Discussion 

We have shown that the space density of 
optically-luminous quasars, M J 4s(1450A)< — 26™9, 
declines moderately, if at all, from z w 2 to z w 4, 
contrasting sharply with the steep decline (factor 
~ 10) found for samples including lower luminos- 
ity quasars (Shaver et al. 1996, 1999, Warren et 
al. 1994, Schmidt et al. 1995, Kennefick et al. 
1995). 

We stress that the quasars in our sample come 
from the bright tail of the optical luminosity func- 
tion. The difference between our result and those 
previously reported is due to the different cosmo- 
logical histories of quasars of different optical lu- 
minosity, the space density of the highly- luminous 
ones having peaked at earlier epochs than that of 
the less luminous ones. 

This abundance of very bright quasars at z w 4 
has already been noted by Fan et al. (2001b), who 
showed that their luminosity function at z ~ 4 
is flatter than that at lower redshifts (note that 
this is not readily- apparent in their Fig. 4, where 
the z ~ 2 2dF-luminosity function is not reported 
correctly) . Our z = 4 cumulative luminosity func- 
tion is broadly consistent with that of Fan et al. 
(2001b) for SDSS, but extends it by about one 
decade in optical luminosity. 

The relatively high space density of optically- 
luminous quasars at high redshift, and the con- 
sequent existence of super-massive black holes at 
very early stages of the evolution of the uni- 
verse, yield non-trivial constraints on models for 
structure formation (e.g.: Haehnelt & Rees, 1993; 
Granato et al. 2001; Haiman & Locb 2001). 

6. Conclusions 

We sought z > 3.8 quasars identified with 
FIRST radio sources Siaghz > 1 niJy (irre- 
spective of radio spectral index), i.e. 30 times 
fainter than previous searches, and built a sam- 
ple of 13 optically-bright radio-loud quasars 
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(M AB (1450A)< -26 m 9, logPi^G/f^WHz" 1 ) > 
25.7) with 3.8 < z < 4.5. Our conclusions are: 

• The space density of quasars meeting the 
above luminosity criteria is (1.0 ± 0.3) 
Gpc- 3 . 

• From a comparison sample at z m 2 we find 
that the space density of radio-loud quasars 
declines between z ~ 2 and z ~ 4 by a factor 
1.9±0.7. 

• We obtained an improved measurement of 
the radio-loud fraction of quasars in the 
luminosity range logPi. 4 (WHz" 1 ) > 25.7: 
Grl = (13.4 ± 3)%, nearly independent of 
optical luminosity, in the range — 25 m > 
M AB (1450A) > -29 m . 

• Using the above value of Grl, we calcu- 
late the space density of optically-selected 
quasars with M AB (1450A)< -26 m 9, to be 
(7.4 ± 2.6) Gpc" 3 . Using the 2dF quasar lu- 
minosity function, we find a ratio p Z 2 / p Z 4 = 
1.8 ± 0.8, consistent with that for the radio- 
loud quasars, but significantly different from 
the generally quoted ratio ~ 10 found for 
samples including lower-luminosity quasars. 

• The moderate decline of quasar space den- 
sity from z = 2 to z = 4 may halt altogether 
at the very bright end of the optical luminos- 
ity function, as, to some extent, is apparent 
in results from SDSS (Fan et al. 2001b). 



Appendix 
putation 



Absolute magnitude com- 



We obtained the continuum absolute magni- 
tudes of the quasars, M AB (1450A) using: 

M AB (1450A) = E - 43.264 - 51og(d L (z)) 

-2.51og/ ^ oga (A)T(A)dA - 2.51og(l + z) 

+2 .51og f g ?(V(i+ ? )) T / A ) dA 

E is the APS POSS-I magnitude corrected for 
Galactic extinction. AB apparent magnitudes are 
defined by m A e(A) = — 2.51og5 , j / (A) — 48.60, where 
S v is measured in [erg cm~ 2 s _1 Hz -1 ] (Okc 
1974). S^ oga in [erg cm" 2 s" 1 A" 1 ] was taken 



- 1.5 

< 

5 







i 





Fig. 5. — The top panel shows the FBQS com- 
posite quasar spectrum (Brotherton ct al. 2001), 
overplotted with a power law with spectral index 
a v = —0.5 (dotted curve). The bottom panel 
shows the quasar template (solid curve) in the ob- 
server's frame for z = 4.2, after correction for the 
Madau (1995) transmission. The POSS-I E re- 
sponse is shown as a dashed curve. 

from Fukugita et al. (1996). is based on the 
FBQS composite by Brotherton et al. (2001) and 
on the models of Madau (1995) for intervening ab- 
sorption. First, an unabsorbed FBQS composite 
was created using the Madau model for z = 2.5. 
This was then corrected for intervening absorp- 
tion at each redshift using linear interpolations of 
Madau's transmission curves for z = 3.5 and z 
= 4.5 (see Fig. 5 for an example). T(A) is the 
response function of the POSS-I E filter + plate 
(Evans 1989) and dh is the luminosity distance in 
Mpc. The E magnitude of Vega, P(Vega) = 0^05, 
was obtained from the relation E — R vs V — R 
by Humphreys et al. (1991), using U(Vcga) and 
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Fig. 6. — The if -correction between apparent 
magnitude m J 4s[1450(l+z)] and apparent E mag- 
nitude (POSS-I plate+filter response). The solid 
line corresponds to Lya+NV equivalent width 70 
A, the dotted line to 15 A. 

i?(Vega) from Fukugita et al. (1995). 

The combined Lya+NV rest-frame equivalent 
width (EW) of the FBQS composite is 70 A at 
z = 4.2, close to the mean, 75 A (rms <~ 30 A), 
for our z = 4 sample. The Ly a line of our quasar 
0918+0636 is heavily extinguished, EW(Ly a + 
NV)< 15 A. For this quasar the emission-line con- 
tour of the template was scaled to EW = 15A. 

The relation between toab [1450(1 + z)\ and E 
is shown as a function of redshift for 3.8 < z < 4.5 
in Fig. 6 for the two values of EW. E magnitude 
decreases with increasing strength of the emission 
line, particularly at z ~ 4.2, when Lya+NV fall 
in the middle of the E band. 
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